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EXECUTIVE  SUMMARY 


OYSryllfil 


Svvaral  rscant  studiM  hav«  shoMn  that  bom  of  tho 
fundaaontal  asauaptiona  uaad  to  darivo  tha  raquiraaanta 
coaputation  foraulaa  currantly  uaad  in  tha  0062  Econoaic  Ordar 
Quantity  (EOG)  Buy  Coaputation  Syataa  provida  a  poor 
approxiaation  to  tha  actual  char actor i at ica  of  AFLC  EOB  itaaa. 
Thia  papar  auaaarizaa  tha  raaulta  of  a  atudy  to  davalop  iaprovad 
daacriptiona  of  EOQ  daaand  and  load  tiaa  char ac tar i at ica,  and  to 
dariva  and  avaluata  altarnata  raquiraaanta  calculationa  baaad 
upon  thia  naa  knowladga. 

APPROACH 

Hiatorical  data  »«aa  collactad  and  analyzad  diacribing  0062 
itaa  daaand  hiatoriaa,  raquiaition  aizaa  and  prioritiaa,  load 
tiaa  variablility,  and  aaaociatad  aircraft  flying  hour  prograaa. 
Ha  than  davalopad  fiva  now  raquiraaanta  coaputationa  for 
dotailad  coat-affactivanaaa  coapariaona  with  tha  currant  0062 
rulaa.  Each  of  thaaa  now  coaputationa  taka  difforont  approachaa 
in  raaolving  tha  aavaral  atatiatical  and  analytic  aodaling 
iaauaa  aaaociatad  with  tha  davol opaant  of  a  practical 
raquiraaanta  aathodology.  Tha  Invantory  Syatoaa  Siaulator 
(INSSIM)  Maa  than  uaad  to  aiaulata  how  aach  of  tha  propoaad 
aathoda  would  hava  parforaad  had  tha  now  aathod  boon  uaad  to 
oanaga  0062  itaaa  during  tha  CY73-79  intorval.  It  waa  found 
that  four  of  tha  fiva  propoaad  policiaa  providod  aignificant 
coat-affactivanaaa  i apr ovaaanta  ovar  tha  currant  0062  aathoda. 


magi  EiNfiiiss 

1.  0062  daaand  pattarna  ara  auch  aora  arratic  than  ia 
aaauaad  in  currant  0062  raquiraaanta  foraulaa.  Thia  aaana  that 
aignificantly  highar  aafaty  atocka  ara  naadad  to  prcivida  high 
lavala  of  aupply  aupport  than  ara  paraittad  in  tha  currant 
coaputation. 


2.  Load  tiaa  variablility  aignificantly  incraaaaa 
uncartainty  aaaociatad  with  itaa  laad  tiaa  daaanda.  Howa> 
tha  currant  0062  rulaa  aaauaa  that  laad  tiaa  variability 
naglagibla. 


3.  Foracaat  arrora  for  F— 104  and  F— 3  aircraft  itaaa  diffar 
aignificantly  froa  thoaa  aaaociatad  with  othar  itaaa  in  tha 
INSSIM  Oata  Bank.  Ha  racoaaand  that  F>104  and  F-5  itaaa  ba 
aaparataly  analyzad  in  futura  INSSIM  aiaulation  atudiaa. 


4.  Raquisition  count  data  containad  in  tha  D062  ayataa 
appaars  to  ba  vary  unraliabla,  and  aa  balittva  that  any 
calculation  basad  on  this  data  — including  tha  currant  0062 
aafaty  laval  calculation —  is  of  littla  valua.  Wa  baliava  that 
iaaadiata  actions  should  ba  takan  to  aithar  (a)  corract  tha  data 
systaa  problaas  associatad  with  raquisitcxi  siza  statistics,  or 
(b>  iaplaaant  altarnata  raquiraaants  foraulas  that  do  not 
raquira  itaa  raquisition  siza  statistics. 

5.  Laad  tiaa  variidiility  has  aajor  iapacts  upon 
raquiraaants  for  safaty  stocks.  Howavar,  accurata  data 
dascribing  laad  tiaa  variability  is  axtraaaly  difficult  to 
obtain  froa  currant  AF  data  systaas.  Additional  work  to  iaprova 
capabilitias  for  laad  tiaa  foracasting  is  graatly  naadad. 

6.  Fcxir  naw  coaputation  policias  wara  idantifiad  which 
appaar  to  ba  significantly  aora  cost-af f activa  than  tha  currant 
0062  coaputations.  Howavar,  bacausa  of  tha  data  problaas 
dascribad  in  (4)  and  <S)  abova,  only  ona  of  thasa  tachniquas  is 
a  candidata  for  iaaadiata  iaptcNMaitaticm. 


REC0WMEN0AT10W3 


1.  Ha  racoaaand  that  Policy  Coda  80,  dascribad  in  Appandix 
B  of  this  raport,  ba  considarad  for  iaaadiata  iaplaaantation  in 
tha  0062  systaa.  This  coaputation  policy  appaars  to  ba 
significantly  aora  cost  *-af  fact  iva  than  tha  currant  0062 
foraulas,  and  it  can  ba  iaplaaantad  using  currant  0062  data. 
Furthar,  only  a  faw  linas  of  coaputar  coda  naad  to  ba  changad  to 
adopt  thasa  calculations.  Howovar,  bafora  iaplaaantation 
afforts  ara  startad,  wa  racoaaand  that  EOQSin  ba  usad  to 
avaluata  tha  transiant  affacts  of  iaplaaanting  this  policy. 
Analytic  calculations  indicata  that  tha  naw  policy  usas  a 
significantly  diffarant  stratagy  in  satting  safaty  stocks  than 
tha  currant  0062  foraulas.  This  naans  that,  unlass  tha  phasa-in 
is  carafully  aanagad,  a  larga  surga  of  buying  activity  could 
occur  whan  tha  naw  foraulas  ara  first  adoptad.  Tha  EOQSIH 
projactions  should  quantify  tha  axtant  of  this  problan,  and 
should  provida  a  basis  for  tasting  altarnata  iaplaaantation 
stratagias. 

2.  Estiaatas  of  0062  itaa  laad  tiaas  hava  a  draaatic 
affact  upon  EOQ  raquiraaants  coaputations,  yat  it  is  vary 
difficult  to  avaluata  tha  accuracy  or  variabilility  of  laad  tiaa 
foracasts  using  currant  AF  data  systaas.  To  provida  a  basis  for 
futura  iaprovaaants  in  raquiraaants  coaputations,  wa  racoaaand 
that  routina  aathods  ba  davalopad  for  aaasuring  tha  variability 
of  itaa  laad  tiaa  and  tha  accuracy  of  laad  tiaa  foracasts,  and 
for  raporting  this  inforaation  within  tha  0062  systaa. 
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ON/ERVIEM 

At  pmcnt,  AFLC  EccNioaic  Ordar  Quantity  (ECXl)  inventory 
•anagMMKit  policies  are  based  upon  eany  of  the  coeeonly  accepted 
assueptions  described  in  inventory  eanageeent  texts,  and 
sieulation  studies  have  sho»«n  that  these  methods  are 
significantly  more  cost-effective  than  previously  used  policies. 
However,  several  recent  studies  have  shown  that  several  of  these 
fundamental  assumptions  provide  poor  d^proxi mat ions  to  actual 
AFLC  item  characteristics.  In  this  paper,  we  summarize  the 
results  of  a  study  to  develop  improved  understanding  of  the 
characteristics  of  AFLC  EOQ  items,  and  to  develop  and  evaluate 
alternate  inventory  management  policies  which  exploit  this 
improved  knowledge.  Several  policies  were  identified  which 
appear  to  be  more  cost-effective  than  current  AFLC  D062 
requirements  computations. 


STUDY  OBJECTIVES: 

This  study  had  three  major  objectives: 

1.  To  develop  improved  understandings  and 
characteristics  of  Economic  Order  Quantity  demand 
processes 

2.  To  develop  alternate  requirements  computation 
techniques  using  a  more  precise  description  of  EOQ 
item  characteristics  than  is  employed  in  current 


methods 


3 


To  coaparo  tho  cost.-offoctivonomm  of  tho  noMly 
dovolopod  procoduros  with  curront  D062  policia 


APPROACH 

To  accoaplish  tha  4^ova  objactivas,  i«a  bagan  by 
parforaing  a  nuabar  of  statistical  analysas  of  actual  EOQ  daaand 
historias  containad  in  tha  INSSIH  data  bank.  This  data  bank 
contains  actual  daaand  historias  for  approxiaataly  40,000  0062 
itaas  aanagad  by  tha  Sacraaanto,  Oklahoaa  .City,  and  Warnar  Robins 
Air  Logistics  Cantars  (ALC) .  Thasa  historias  covar  tha  CY71*-CY79 
intarval,  a  total  of  38  quartars  of  daaand  data.  A  aajor 
objactiva  of  thasa  statistical  studias  Mas  to  davalop  an 
aapirical  aodal  of  tha  distribution  of  foracast  arrors  for  thasa 
itaas.  Wa  also  collactad  statistical  data  dascribing  itaa 
raquisition  sizas,  raquisition  prioritias,  and  itaa  laad  tiaa 
variability.  This  lattar  inforaation  Mas  priaarily  obtainad  froa 
pravious  AFLC  study  af forts. 

Araad  Mith  a  aora  pracisa  dascription  of  tha 
char actari sties  of  EOQ  itaas,  aa  davalopad  savaral  altarnata 
analytical  aodals  of  EOQ  daaand  procassas,  and  ua  than  darivad 
tha  asscxriatad  "optiaua"  raquiraaants  coaputation  policias.  Each 
of  thasa  aodals  provida  a  aora  datailad  dascription  of  actual  EOQ 
itaa  char actari sties  t^ ^n  ara  jtilizad  in  tha  currant  0062 
raquiraaants  coaputation.  r,jMavar,  aach  altarnata  aodal  takas  a 
diffarant  approach  to  tha  ceabination  of  statistical  astiaation 
and  aathaaatical  aodal ing  issuas  involvad  in  tha  davalopaant  of  a 
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practical  requiraaenta  coaputation  prcxrass. 

In  tha  final  staga  of  thia  atudy,  wa  parforaad  data! lad 
aiaulation  axpariaanta  to  coapara  tha  ralativa  coat— af f activanaaa 
of  tha  naaly  davalopad  policiaa  with  tha  currant  D062 
raquiraaanta  coaputationa.  In  parforaing  thaaa  aiaulation 
atudiaa  wa  ccjnatructad  four  aaparata  itaa  aaaplaa  of 
approxiaataly  fiva  hundrad  itaaa  aach  froa  inforaation  in  tha 
INSSIM  Data  Bank.  Wa  than  uaad  tha  Invantory  SyataM  Siaulator 
(INSSIM)  to  avaluata  aach  of  tha  propoaad  policiaa.  Basically, 

•  p 

our  aiaulation  axpariaanta  avaluated  how  each  of  tha  proposed 
policies  would  have  parforaad  in  aanaging  aach  of  thaaa  sets  of 
itaaa  during  the  CY73-79  interval.  Froa  these  coaputar  runs,  a 
lar-a  nuabar  of  coat-aff activanaaa  curves  ware  constructed  and 
analyzed.  As  a  result,  we  identified  four  new  policies  which 
appear  to  be  superior  to  tha  currant  0062  raqui r eaant s 
coaputation  rules.  One  of  these  appears  to  be  particularly 
suitable  for  iaaediata  iaplaaentation. 

Finally,  wa  docuaentad  the  results  of  this  study  in  a 
series  of  eight  detailed  working  papers.  These  reports  are 
listed  in  Appendix  A.  In  this  paper,  wa  .  suaaariza  tha  aajor 
results  of  tha  study  effort. 

With  this  background,  let  us  now  discuss  the  aajor 
findings  and  conclusions  obtained  froa  this  study.  For 
convanianca,  wa  have  organized  our  discussion  into  three  aajor 
sections!  these  correspond  to  tha  Statistical  Analysis,  Analytic 
Modeling,  and  Siaulation  aspects  of  tha  study  effort. 
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STATISTICAL  RESULTS 

As  noted  above,  we  began  our  e-ffort  by  performing  a 
number  of  statistical  studies  to  develop  improved  understanding 
of  the  characteristics  of  EOQ  items  and  of  the  sources  of 
forecasting  errors  associated  with  these  items.  We  began  our 
statistical  efforts  by  collecting  actual  and  predicted  flying 
hour  data  for  each  of  the  twenty-three  aircraft  represented  in 
the  INSSin  data  bank.  We  then  computed  an  "Accuracy  Ratio",  the 
ratio  of  total  flying  hours  predicted  for  a  given  twelve  month 
interval  to  by  the  associated  actual  flying  program  for  that 
interval.  For  example,  as  of  1  January  1973,  a  total  of  230,800 
flying  hours  were  forecast  for  B— 52  aircraft  for  the  interval 
January  through  December  1973.  In  fact,  a  total  of  203,703  hours 
were  flown  by  the  BS2  during  this  interval.  Hence,  the  accuracy 
ratio  associated  with  this  forecast  equaled  (230,800) / (203,703)  * 
1.133.  That  is,  the  January  forecast  exceeded  the  eventually 
observed  program  by  13. 3X. 


PA6E  S 


Similar  calculations  Mmrm  pmrformmd  -for  mach  -forecast 
made  during  the  CY73-CY79  interval,  and  for  mach  of  thm  aircraft 
represented  in  the  INSSIM  data  bank.  The  results  of  these 
calculations  are  presented  in  Figure  1.  We  then  performed 
Analysis  of  Variance  (ANOVA)  tests  to  determine  if  there  was  a 
statistically  significant  difference  among  aircraft  or  among  time 
periods  with  respect  to  flying  program  forecast  accuracy.  After 
deleting  observations  asociated  with  very  low  levels  of  flying 
activity,  we  found  no  statistically  significant  difference  in 
flying  program  forecast  accuracy  performance  among  weapons  or 
across  time  frames.  Consequently,  it  appears  reasonable  to  pool 
these  estimates  to  describe  flying  program  forecast  accuracy. 
This  pooled  distribution  is  presented  in  Figure  2.  As  shown  in 
the  figure,  flying  program  accuracy  ratios  range  from  76.3 
percent  to  181.2  percent  of  the  forecasted  values,  and  forecasts 
averaged  12  percent  above  the  flying  program  that  was  eventually 
observed.  Hence,  ii«e  conclude  that  errors  in  flying  program 
forecasts  are  a  potentially  significant  source  of  forecasting 
errors  in  EOQ  requirements  projections. 

Next,  we  performed  a  number  of  correlation  and  Analysis 
of  Variance  studies  to  develop  an  improved  description  of  the 
relationship  between  forecast  EOQ  usage  rates  and  the  actual 
demands  for  these  items.  In  our  initial  efforts,  we  asked  if 
there  were  any  significant  differences  in  the  magnitude  of  demand 
forecast  errors  across  aircraft,  across  item  demand  rate 
categories, across  Air  Logistics  Centers,  or  across  time  periods. 


Figure  I.  Forecast  Accuracy  Ratios  for  CY73-Cy79 
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During  these  periods^  actual  flying  hours  were  so  low  that 

tre  accuracy  ratios  are  not  meaningful.  To  prevent  these 
points  from  biasing  our  results  upward,  the  true  values 
were  replaced  by  1.000. 
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Figure  2. 

Frequency  Distribution  of 
Forecast  Ratios,  CY  73-79. 
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W«  found  that,  yos,  significant  diffarancas  axist.  In 
particular,  Na  found  that  foracaating  arrors  aaaociatad  with 
F-104  and  F-5  aircraft  and  foracasta  aaaociatad  with  tha  CY75-76 
tiaa  intarval  wara  significantly  diffarant  in  aagnituda  froa 
thosa  aaaociatad  with  tha  othar  21  aircraft  in  tha  INSSIM  data 
bank  or  with  othar  tiaa  parioda.  Subaaquantly,  wa  found  that  tha 
F-104  was  phasad  out  of  tha  USAF  invantory  during  tha  CY7S-76 
intarval,  but  was  latar  raintroducad  as  a  aajor  Foraign  Military 
Salas  aircraft.  Vary  significant  foracaating  arrors  wara 
associatad  with  this  phasa— out/phasa— in  procass.  On  tha  othar 
hand,  wa  found  that  tha  F-S  is  tha  only  aircraft  in  tha  INSSIM 
data  bank  that  has  a  continuously  incraasing  prograa. 
Unfortunataly,  thara  ara  vary  faw  F-5  itaas  in  tha  INSSIM  data 
bank,  so  tha  construction  of  a  singla  data  sat  consisting  of  only 
F-5  itaas  is  not  a  practical  analysis  altarnativa.  Consaquantly, 
it  is  racoaaandad  that  in  futura  siaulation  studias,  F-104  and 
F-5  itaas  should  ba  saparataly  analyzad  froa  tha  othar  INSSIM 
aircraft.  Furthar,  it  is  racoaaandad  that  C-5,  F-15,  F-16,  or 
siailar  0062  itaas  ba  addad  to  tha  INSSIM  Data  Bank  to  parait 
futura  analysis  of  itaas  with  ascanding  prograas. 

Nith  tha  abova  inforaation,  wa  again  rapaatad  tha 
corralation  and  ANOVA  studias  with  all  F-104  and  F— 5  itaas 
dalatad.  Bacausa  of  tha  larga  nuabar  of  itaas  in  our  data  sat, 
statistically  significant  diffarancas  wara  still  obsarvad. 
Howavar,  tha  diffarancas  across  aircraft  and  across  tiaa  pariods 
wara  of  saall  aagnituda  and  appaar  to  ba  of  littla  practical 


PAGE  9 


impact.  Consequently,  we  believe  that  statistical  estimates  of 
forecast  errors  may  be  pooled  across  aircraft  and  time  periods 
for  the  remaining  IN5SIM  items.  However,  major  differences  in 
error  performance  still  exist  among  high  and  low  activity  item 
classes.  Consequently,  separate  empirical  models  were  developed 
for  high  and  low  activity  item  classes. 

Next,  we  performed  several  studies  to  identify  the  shape 
of  EOQ  demand  rates  forecast  error  distributions.  We  found  that, 
in  general,  the  demand  processes  for  EOQ  items  are  much  more 
erratic  than  is  assumed  in  the  current  D062  formulas.  At 
present,  EOQ  demand  rates  are  forecast  by  first  estimating  the 
average  usage  rate  over  the  past  two  years,  and  then  adjusting 
this  value  for  anticipated  changes  in  flying  program  activity.  It 
is  then  assumed  that  actual  demand  will  be  normally  distributed 
about  this  forecasted  rate.  That  is, in  D062  it  is  assumed  that 
the  actual  demand  will  be  distributed  about  the  forecasted  value 
according  to  the  familiar  "bel 1 -shaped"  curve.  In  contrast,  we 
found  that  the  distribution  of  actual  quarterly  demand  about  the 
forecasted  value  is  a  highly  skewed  distribution  with  significant 
levels  of  probability  in  the  right-hand  tail  of  the  distribution. 
This  means  that  significantly  higher  levels  of  safety  stock  are 
required  to  provide  a  given  fill  rate  than  would  be  required  if 
the  normal  probability  model  was  in  fact  true. 

For  example.  Figure  3  compares  the  normal  distribution 
with  observed  distributions  of  EOQ  demand  forecast  errors.  The 
four  distributions  shown  describe  the  relationship  of  actual 


Figure  3.  An  Observed  Distribution  o£  Period  Errors 
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demands  to  the  corresponding  forecasted  values  for  forecasts  of 
one,  tMO,  three,  and  four  periods  into  the  future.  The 
standardized  error  Zi  is  defined  as 

Actual  Demand  _  Forecasted  Demand 
in  Period  i  for  period  i 

Zi  = 

MAD 

where  MAD  denotes  the  Mean  Absolute  Deviation  of  past  demand, 
i.e.  the  average  magnitude  of  past  usage  about  the  average  value. 
Thus  Z3  denotes  the  distribution  of  actual  demnd  in  a  quarter 
relative  to  the  usage  rate  that  was  forecast  three  quarters 
earlier.  As  shotet  in  the  Figure,  if  demand  is  normally 
distributed  there  is  a  SO  X  chance  that  actual  demand  will  be 
less  than  or  equal  to  the  forecasted  value,  and  it  is  almost 
certain  that  actual  demand  will  be  lass  than  four  MADs  above  the 
forecasted  value.  In  contrast,  we  observed  that  for  low  demand 
items  there  is  almost  an  80  X  chance  that  actual  demands  will  be 
less  than  forecasted  value,  and  there  is  at  least  a  5  X  chance 
that  demands  for  these  items  will  exceed  four  MADs  above  the 
forecast.  Almost  identical  curves  are  obtained  for  all  time 
periods  and  across  all  Air  Logistics  Center  item  samples.  Fcsr 
high  activity  items,  the  observed  error  distributions  are  closer 
to  the  normal  curve,  but  the  error  distributions  are  still 
distinctly  skewed.  As  noted  above, the  observed  error  distribution 


PAGE  12 


iapliM  that  if  high  sarvica  lavala  ara  to  ba  obtainad, 
significantly  highar  safaty  stocks  ara  raquirad  than  is  naadad 
»«han  foracast  arrors  ara  norally  distributad. 

Tha  above  data  describes  the  distribution  of  deaand  in  a 
single  period.  Wa  also  studied  tha  distribution  of  total  net 
deaand  observed  in  a  lead  tiaa  one,  two,  .  .  .  ,  or  eight 
quarters  relative  to  tha  D062  foracast  usage  for  tha  combined 
iterval.  Wa  originally  anticipated  that  tha  resulting 
distribution  mould  ba  described  by.  one  of  tha  standard 
probability  models  found  in  statistics  and  inventory  management 
texts.  Homavar,  tat  found  that  tha  observed  lead  time  demand 
distribution  does  not  fit  any  of  these  standard  models. 
Consequently,  me  developed  reasonable  analytic  approximations  for 
the  observed  distribution  of  demand  in  a  lead  time.  These 
approximations  use  exponential  functions  to  approximate  the 
observed  lead  time  demand  distribution.  We  refer  to  the 
resulting  approximation  as  the  exponential '-constant,  or  EXP-CON 
model,  since  this  model  is  simply  an  exponential  approximation  to 
the  distribution  of  demand  in  a  lead  time  that  is  a  knomn 
constant.  A  comparison  of  the  normal  and  EXP-CON  models  with  the 
observed  distribution  of  demand  is  presented  in  Figure  4. 

Analysis  of  lead  time  variability  data  was  another 
important  research  area.  The  above  studies  quantify  the 
variability  of  a  demand  in  a  given  number  of  periods.  When 
replenishment  lead  times  are  also  variable,  the  distribution  of 
demand  in  a  lead  tine  is  even  more  variable  than  is  illustrated 


Figure  4.  Comparison  of  the  Normal  and  Exponential 
with  Observed  Demand  Distribution. 
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in  Figures  3  and  4. 

We  -found  that  data  describing  lead  tiee  variability  is 
extreeely  difficult  to  obtain.  The  best  source  of  infornation 
Mhich  Me  could  locate  is  contained  in  a  1981  study  by  Dr.  Jack 
Hayya.  Hayya  studied  the  variability  of  replenish«»ent  lead  tiees 
for  62  high  activity  D062  itees  and  attempted  to  fit  several 
probability  distributions  to  the  observed  lead  tiee  data.  Based 
on  Hayya’s  results.  Me  believe  that  the  Baema  probability 
distribution  provides  the  most  reasonable  single  model  of 
replenishment  lead  time  variability. 

The  Gamma  distribution  is  illustrated  in  Figure  5  using  a 
coefficient  of  variation  of  .333.  This  is  the  median  coefficient 
of  variation  value  for  the  items  studied  by  Hayya.  That  is,  half 
of  the  items  studied  by  Hayya  had  less  variability  than  indicated 
in  this  figure,  Mhile  approximately  half  had  more  variability. 
We  used  this  model  of  lead  time  variability  as  the  basis  for  our 
simulation  studies.  We  believe  that  this  is  an  adequate  model 
for  use  in  this  study.  Hoeever,  more  research  is  required  if 
accurate  lead  time  variability  estimates  are  to  be  developed  for 
all  D062  items. 

Finally,  if  we  assume  that  demands  in  a  given  lead  time 
are  described  by  the  exponential  approximation  illustrated  in 
Figure  4,  and  if  we  also  assume  that  replenishment  lead  time  is 
independent  of  demand  and  described  by  a  Gamma  distribution  such 
as  that  shown  in  Figure  S,  then  analytical  methods  may  be  used  to 
compute  the  probability  distribution  of  demand  in  a  lead  time. 


Figure  5.  Gamma  Density  Function  with  Mean 


PAGE  16 


We  refer  to  the  resulting  distribution  as  the  Exponential-Gaeea 
lead  time  demand  model,  or  EXP-G^  -for  short. 

The  relationships  among  the  Normal,  EXP-CON,  and  EXP-GAM 
lead  time  demand  models  are  illustrated  in  Figure  6.  As  shown  on 
the  -figure,  the  EXP-GAM  lead  time  demand  distribution  is  even 
more  erratic  than  the  EXP— CON  a»del.  Thus,  even  higher  safety 
stocks  are  required  to  provide  high  service  levels  in  this  latter 


Standardized  Demand  In  Leadtime 
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ANALYTICAL  MODELING 

During  th»  analytical  aodaling  phasa  o-f  thie  project  we 
sought  to  develop  new  EOQ  reorder  point  and  order  quantity 
calculations  based  upon  a  more  accurate  description  o-f  EOQ  itee 
characteristics  than  is  employed  in  current  D062  formulas.  In 
this  effort,  we  developed  several  alternate  models  of  EOQ  demand 
processes,  and  we  then  derived  the  associated  “optimum" 
requirements  computation  policies.  As  noted  above,  each  of  these 
models  provide  a  more  detailed  description  of  actual  EOQ  item 
characteristics  than  current  formulas,  but  each  of  the  alternate 
models  resolve  important  statistical  estimation  and  analytic 
modeling  issues  in  different  ways. 


Basically,  three  categories  of  inventory  management  policies 
developed  for  further  testing.  These  aret 


1.  The  current  D062  policy  and  modifications  of  this  policy 
to  account  for  demand  variability  and  for  leadtime  variability. 

2.  Policies  derived  from  knowledge  of  the  distribution  of 
units  demanded  per  requisition  and  of  the  distribution  of  item 
lead  times. 


3.  Policies  derived  from  the  observed  distribution  of 
demand  in  a  given  lead  time  obtained  from  aggregate  error 
statistics. 


f 

*■ 
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Frttm  th«s«  three  basic  categories,  six  policies  Mere 
selected  for  detailed  testing.  These  basic  policies  are 
described  in  Table  I. 

As  shoMn  in  Table  I,  each  policy  Mas  assigned  a  "Policy 
Code"  to  sieplify  discussions  of  the  basic  techniques.  Policy 
Code  10  refers  to  the  current  0062  formulas.  In  studying 
individual  item  demand  patterns.  Me  observed  that  many  item 
demand  histories  contain  "spikes",  i.e.  very  large  erratic 
demands  that  are  inconsistent  Mith  demands  either  before  or  after 
the  spike.  Hence,  Policy  Code  20  is  identical  to  Policy  lO, 
except  that  "spikes"  are  deleted  froq  demand  rate  estimates  to 
prevent  an  undesireable  upMard  bias.  Similarly,  Policy  Code  60 
is  also  identical  to  the  Policy  lO  formulas,  Mith  the  exceptici* 
that  the  standard  deviation  estimate  is  modified  to  account  for 
lead  time  variability. 

In  the  second  category  of  models.  Me  began  our  analyis  by 
considering  a  situation  in  Mhich  (a)  customers  arrive  according 
to  a  Poisson  process  Mhose  mean  is  proportional  to  flying  program 
activity,  (b)  requisiton  sizes  obey  a  logarithmic  distribution, 
<c>  lead  times  are  independent  random  variables  Mhich  satisfy  a 
gamma  probability  distribution,  <d)  all  requisitions  are  assumed 
to  be  of  equal  priority,  and  (e>  a  continuous  revisM  (Q,R)  order 
quantity,  reorder  point  control  system  is  being  employed.  With 
these  assumptions,  the  distribution  of  demand  in  a  lead  time  is 
described  by  the  Logarithmic-Poisson-Gamma  (LP6>  distribution. 
This  is  a  noM  probability  distribution  Mhich  Mas  derived 
•specially  for  this  study. 
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Table  1 

Inventory  Management  Policy  Codes 

Inventory  Management  Policy 
Current  D062  Formulas 

Current  D062  Formulas,  with  outliers  excluded  from  demand  and  vari* 
ance  estimates 

Current  D062  Formulas,  with  adjustments  to  standard  deviation  of 
lead  time  demand  to  account  for  lead  time  variability 

Scaled  Negative  Binomial  reorder  point  calculations 
Constant  Leadtime  Exponential  Forecast  &ror  model 
Exponential-<iamma  Forecast  Error  Model 
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Unfortunately,  tha  LP6  distribution  is  quit*  coaplsx,  and 
significant  conputaticwial  difficulties  are  encountered  in 
evaluating  percentage  points  of  this  distribution  for  other  than 
vary  Iom  values  of  expected  lead  tiae  deeand.  At  the  beginning 
of  this  study,  we  felt  that  the  LP6  distribution  provided  the 
most  feasible  detailed  analytic  description  of  AFLC  EOQ  deeands. 
Because  there  are  over  500,000  itees  managed  by  the  0062  system, 
we  considered  the  development  of  an  efficient  computational 
approach  for  the  LP8  distribution  to  be  a  primary  research 
objective.  Consequently,  significant  efforts  ware  devoted  to 
develop  computationally  efficient  approximations  to  the  LP6 
amdel.  We  found  that  a  Scaled  Negative  Binomial  probability 
distribution  provides  the  required  approximation.  Details  of  the 
Scaled  Negative  Binaomial  probability  model  are  presented  in 
Working  Paper  81-04.  On  the  other  hand,  computer  programs  for 
computing  percentage  ponts  for  the  LP6  and  Scaled  Negative 
Binomial  distributions  are  presented  in  Working  Paper  81-05. 
Policy  Code  70  described  in  Table  Z  corresponds  to  the  use  of  the 
Scaled  Negative  Binomial  model  for  requirements  computations. 

We  know  from  sensitivity  analyses  of  the  LP6  distribution 
that  the  distributon  of  individual  requisition  sizes  has  an 
important  influence  on  the  shape  of  the  lead  tiae  demand 
distribution.  Since  the  LP8  model  has  four  parameters,  it  is 
possible  to  more  closely  match  the  LPO  model  to  the  specific 
statistical  estimates  associated  with  a  given  item  than  is 


possibls  using  othsr  ssthcxis.  Un-fortunatsly*  this  Matching  is 
not  a  frmm  good.  Bacausa  of  tha  liaitad  data  avai labia  for  any 
singla  itam,  singla  itaa  astiaatas  ara  subjact  to  significant 
statistical  astiaation  arror.  Ona  aathod  to  raduca  thasa 
statistical  problaas  is  to  usa  aggragata  itaa  statistics.  In  tha 
naxt  paragraph,  wa  discuss  policias  basad  on  this  lattar 
approach. 

Tha  final  catagory  of  analytical  acjdals  Mhich  aa  studiad 
utiliza  daaand  foracast  arror  statistics  aggragatad  ovar  larga 
populations  of  itaas  to  dascriba  tha  distribution  of  daaand  in  a 
givan  laad  tiaa.  This  is  tha  Exponantial -Constant,  or  EXP-CON, 
aodal  discussad  abova.  Although  this  aodal  dascribas  aggragata 
error  perforaanca,  it  is  possible  that  no  singla  itaa  is 
accurately  aodaled  by  this  distribution.  Howavar,  by  pooling  tha 
error  statistics  for  a  larga  nuabar  of  itaas,  wa  can  obtain  vary 
precise  statistical  astiaatas  of  this  aggragata  arror 
distribution.  As  shown  in  Table  I,  Policy  Coda  BO  corresponds  to 
tha  EOQ  raquiraaants  coaputations  that  ara  obtained  when  tha  EXP- 
CON  aodal  is  used  to  describe  tha  distribution  of  laad  tiaa 
daaand. 

Tha  final  aodal  selected  for  siaulation  evaluation  is  Policy 
Coda  90.  This  policy  corresponds  to  tha  Exponential -Baaaa,  or 
EXP-6AN,  nodal  of  laad  tiaa  daaand  discussad  above.  Recall  that 
this  nodal  assunas  that  daaand  in  a  givan  laad  tiaa  is  discribad 
by  an  exponential  approx i nation  to  tha  observed  forecast  arror 
distribution,  and  a  gaaaa  distribution  such  as  that  shown  in 


Figura  S  is  ussd  to  dsscribs  Issd  tims  variability.  Nuaaric 
intsgration  tschniquss  arm  than  ussd  to  coaputs  tha  distribution 
of  daaand  in  a  randoa  laad  tiaa. 

Ha  also  attaaptad  to  astand  tha  abova  rasults  to  situations 
Mith  tMo  priority  classas,  but  aa  ancountarad  tao  aajor  problaas. 
First,  aa  aara  unabla  to  locata  any  usaful  sourca  of  individual 
itaa  priority  statistics.  Hanca,  if  a  usaful  analytical  aodal 
aara  davalopad,  aajor  data  systaa  changas  aould  ba  raquirad  to 
support  it.  Sacond,  tha  priority  data  ahich  aa  aara  abla  to 
obtain  indicatas  that  fill  ratas  for  high  priority  raquisitions 
aara  loaar  than  for  routina  priority  raquisitions.  This  suggasts 
that  thara  is  a  coaplicatad  intar dapandanca  aaong  stockaga  policy 
and  tha  priority  of  subaittad  raquisitions.  Thus,  aodals  ahich 
assuaa  indapandimt  priority  straaas  appaar  inconsistant  aith  tha 
availabla  data.  Finally,  avan  ahan  indapandanca  is  assuaad,  tha 
aathaaatical  analysis  laads  to  axprassions  ahich  ara 
coaputationally  impractical  for  usa  in  tha  900, 000+  itam  0062 
systam.  With  sufficiant  affort,  it  may  ba  possibla  to  davalop 
usaful  approximation  aquations  for  this  casai  hoaavar,  aa  fait 
our  afforts  aara  battar  spant  in  improving  tha  LPG  and 
axponantial  approximation  approachas. 
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SIMULATION  RESULTS 

The  Inventory  Systee  Sieulator  CINSSIM)  m«b  used  to  evaluate 
each  of  the  proposed  rules  under  several  different  funding 
levels.  INSSIM  provides  a  detailed  description  of  the  D062 
Econoeic  Order  Quantity  Buy  Coaputation  Systee  and  uses  actual 
Air  Force  demand  histories  to  drive  the  sieulation  process.  For 
this  study,  four  itee  saaples  of  appproxieately  SOO  items  each 
were  selected  from  the  INSSIM  Data  Bank.  The  saaples  SM.H  and 
SM.L  were  selected  froa  Sacraaento  Air  Logistics  Canter  records, 
while  saaples  OC.H  and  OC.L  were  selected  froa  Oklahoma  City  ALC 
records.  The  high  activity  saansles  SM.H  and  OC.H  consisted  of 
items  which  had  net  demands  in  CY71-72  which  exceeded  dSOOO  per 
year,  while  the  saaples  SM.L  and  OC.L  consisted  of  items  with  net 
CY71-72  demands  which  ware  less  than  t  5000  per  year. 

Thirty-eight  quarters  of  history  covering  the  CY71-79 
interval  were  available  for  each  of  these  items.  The  first  eight 
quarters  of  data  were  used  to  initialize  the  forecasting  and 
inventory  management  rules,  while  the  remaining  30  quarters  of 
data  were  used  to  simulate  the  behavior  of  each  of  these  rules. 
Thus,  the  simulation  results  evaluate  how  each  of  these  rules 
would  have  performed  had  they  been  employed  during  the  CY73-79 
interval . 

Each  of  the  policies  discussed  above  were  simulated  under 
several  different  funding  levels,  and  our  results  were  summarized 
in  appropriate  tables  and  cost-effectiveness  plots.  The  detailed 


rasults  of  tho  sioulation  of fort  arm  prosontod  in ^Working  Papor 
81-07.  Briofly,  found  four  policios  Mhoso  sioulation  rosulta 
provided  significant  cost-of focti vonos  i apr ovaosnt s  over  tho 
curront  D062  f  omul  as.  Wo  found  that  Policy  Codes  90,  80,  and  70 
consistently  out-perfore  Policy  Codes  60  and  10  for  the 
requisition  fill  rate,  unit  fill  rate,  average  unit  delay,  and 
average  requisition  delay  eeasures.  Policy  Codes  90  and  80 
perforeed  particularly  eell  for  the  Iom  activity  saeples  and  for 
the  unit-based  measures  of  effectiveness.  On  the  other  hand,  it 
Mas  difficult  to  distingui.sh  a  superior  policy  for  the 
requisition-based  measures  for  the  high  activity  samples. 

Mhen  the  long  supply  versus  buy  dollar  curves  Mere 
considered,  ee  obtained  mixed  results.  Policy  Codes  90,  70,  60, 
and  10  Mere  each  ranked  first  or  tied  for  first  in  at  least  one 
item  sample,  but  Policy  Code  80  Mas  ranked  last  in  all  cases. 
This  is  a  vary  interesting  result,  since  Policy  Code  80  performed 
vary  Mali  Mith  respect  to  each  of  the  fill  rate  and  average  delay 
statistics  reported  above.  Apparently  Policy  Ccxfe  80  achieves 
its  improved  supply-effectiveness  versus  buy  dollar  performance 
by  taking  higher  risks  that  some  of  its  safety  stocks  Mill  later 
be  classified  as  excess.  HoMever,  Policy  Code  80  has  excellent 
support  versus  buy  dollar  curves  on  all  the  other  measures  of 
effectiveness  studied. 
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IHPLICATICMS  FOR  THE  HANAGEHENT  OF  REAL-WORLD  INVENTORIES. 

As  discuBssd  sbovs,  ms  obtsinsd  qusntitstivs  svslustions  for 
ssch  of  ths  six  policiss  in  s  sisulstsd  D062  snvironosnt.  Thsss 
rssults  indicsts  that  Policy  Codss  90,  80,  70,  and  60,  in 
dsscsnding  ordsr  of  prsfarancs,  ars  significantly  aors  cost- 
sffsctivs  than  ths  currsnt  0062  rulss.  Homsvst,  although  ths 
sisulation  providss  a  dstailsd  dsscription  of  ths  D062  systsa, 
thars  ars  ssvsral  iaportant  diffsrsncss  bsttissn  ths  D062 
siaulation  aodsl  and  ths  actual  D062  snvironasnt.  Thsss 
diffsrsncss  aust  ba  considarsd  in  dstqraining  an  appropriats 
policy  for  ths  aanagsasnt  of .ths  actual  D062  systsa. 

First,  in  our  siaulation  aodsl  ths  asan  and  variancs  of 
procursasnt  load  tiass  mots  knosn  Mith  csrtainty,  Mhils  in 
practics  thsss  paraastars  aust  bs  astiaatsd  froa  availabls  data. 
Ws  obssrvsd  that  ssvsral  policiss  significantly  out-parfora  ths 
currsnt  0062  rulss  Mhsn  accurats  Isad  tias  data  is  availabls. 
Unfortunataly,  ths  rsquirsd  Isad  tias  data  is  not  currsntly 
availabls  in  ths  0062  systsa,  and  it  appsars  that  a  si5p)if leant 
data  procsssing  sffort  scxild  bs  rsquirsd  to  routinsly  collact  and 
updats  ths  nssdsd  inforoMtion.  Thus,  iaplaaantation  of  Policiss 
90,  70,  or  60  aust  sithar  <a)  sMsit  ths  dsvslopaant  of  such  a 
systsa  or  (b>  uss  rsgrsssion  or  siailar  sstiaatss  as  an  intaria 
asasurs.  On  ths  othsr  hand.  Policy  10  (ths  currant  0062 
foraulas)  and  Policy  80  do  not  rsquirs  load  tias  variability 
paraastars. 
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Seccmd,  in  C3ur  siaulation  aodel  have  perfect  inforeation 
concerning  itee  requieition  counts.  In  previous  studies  we  have 
found  that  D062  requisition  counts  are  extreowly  unreliable,  and 
we  believe  that  any  forasila  that  uses  0062  requisition  count  data 


— including  the  current  0062  foreulas —  is  basing  Air  Force 
safety  stocks  on  a  randoe  nueber  generator.  This  problee  appears 
particularly  severe  for  the  very  large  nueber  of  low  activity 
itees  managed  by  the  0062  system.  Hence,  we  believe  that  the 
data  processing  rule  of  garbage— in,  garbage-out  wcxild  describe 

the  results  of  implementing  ariy  rule  that  uses  requisition  count 

»  #  ■ 

data  from  the  current  0062  system.  Thus,  rules  70,  60,  and  lO 
ccxild  be  expected  to  perform  much  worse  in  the  "dirty-data” 
environment  of  the  0062  system  than  they  have  performed  in  the 
simulated  system. 


Observe  that  Policy  Code  80  is  the  only  rule  that  avoids 
both  of  the  severe  data  problems  described  above.  Although 
Policy  Ccxfe  90  provided  results  which  were  slightly  superior  to 
Policy  80,  Policy  90  requires  accurate  estimates  of  lead  time 
variability  to  deliver  on  its  promise  of  superior  performance. 


On  the  other  hand.  Policy  80  does  not  require  th 


estimates. 


Policy  80  has  another  clear  advantage.  Hatheoatically,  Policy  80 
is  even  simpler  than  the  current  0062  formulas,  and  only  a  few 
lines  of  code  would  need  to  be  changed  to  implement  this  rule  in 
0062,  in  IN88IM,  in  EOOSIH,  or  in  any  other  data  system  that  uses 
the  current  0062  formulas.  A  flow  chart  documenting  the  logic 
and  calculations  for  Policy  Code  80  is  presented  in  Appendix  B. 
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On  th»  other  hand.  Working  Papar  81-02  docuaents  the  FORTRAN  code 
needed  to  iapleaent  Policy  80  in  IN8SIM,  and  provides  a  logical 
fraeework  -for  the  developaent  of  coeputer  codes  for  other 
effected  systees. 


In  suamary.  Me  believe  that  significant  i eproveeent s  in 
inventory  eanageeent  effectiveness  eay  be  achieved  by  replacing 


the  current  D062  rules  (Policy  10>  by  the  Policy  80  rules  as  soon 
as  possible. 


■f 
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Policy  80  Flowchart 


F>OI  I CY  code;  8 


P(Z  <  k> 


1.000 


exp-con  Dirntribution 


COMPUTE  EOQ 


Q  =  /2AD 
/  IC 


P(k) 


Safety  Factor 

.531  exp (-.463  k> 
for  k  >  O 


BOUND  EOQ 

IF  Q  <  .5D,  set  Q  «  .SO 
IF  Q  >  3D,  set  0  »  3D 


COMPUTE  OPTIMUM  FILL  RATE  F, 

F  -  1  -  0IC/(  7T  D) 

BOUND  FILL  RATE 

IF  F  <  O.,  set  F  -  O 

IF  F  >  .9999,  set  F  -  .9999 


COMPUTE  SAFETY  FACTOR 


A  *  Cost  per  order 
D  *  Annual  Deeand 
Rate  (Units) 

I  *  Holding  (^ost 
1%/%-yr) 

C  B  I tee  Unit  Cost 

*  Shortage  Cost 
per  Unit  Short 

OMAD  ■  Mean  Absolute 
Deviation  (MAO) 
of  Deeand 
per  Quarter 

LT  «  Lead  Tiee  in 
Quarters 


Set  k  *  O 
IF  F  >  .669, 

set  k  -  -Cl/(. 463)3  In  C (1-F) / (. 331 ) 3 


COMPUTE  SAFETY  LEVEL  SL, 


SL  -  k  QMAD  /  LT 
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